Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



On-axis spatially resolved spectroscopy of low redshift quasar 
host galaxies: HE 1503+0228, at z = 0.135 * 

F. Courbin 1 ' 2 , G. Letawe 1 , P. Magain 1 , L. Wisotzki 3 , P. Jablonka 4 , K. Jahnke 5 , B. Kuhlbrodt 5 , D. Alloin 6 , 

G. Meylan 7 , D. Minniti 2 , and I. Burud 7 

1 Institut d'Astrophysique et de Geophysique, Universite de Liege, Allee du 6 Aout, 17, Sart Tilman (Bat. B5C), 
Liege 1, Belgium 

2 Universidad Catolica de Chile, Departamento de Astronomia y Astrofisica, Casilla 306, Santiago 22, Chile 

3 Institut fur Physik, Universitat Potsdam, Am Neuen Palais 10, 14469 Potsdam, Germany 

4 GEPI, Observatoire de Paris, Place Jules Janssen, F-92915 Meudon Cedex, France 

5 Hamburger Sternwarte, Universitaet Hamburg, Gojenbergsweg 112, D-21029 Hamburg, Germany 

6 European Southern Observatory, Casilla 19, Santiago, Chile 

7 Space Telescope Science institute, 3700 San Martin Drive, Baltimore, MD 21218 USA 



Abstract. We present the first result of a comprehensive spectroscopic study of quasar host galaxies. On-axis, 
spatially resolved spectra of low redshift quasars have been obtained with FORS1, mounted on the 8.2m ESO 
Very Large Telescope, Antu. The spectra are spatially deconvolved using a spectroscopic version of the "MCS 
deconvolution algorithm" . The algorithm decomposes two dimensional spectra into the individual spectra of the 
central point-like nucleus and of its host galaxy. Applied to HE 1503+0228 at z = 0.135 (M_g = — 23.0), it provides 
us with the spectrum of the host galaxy between 3600A and 8500A (rest-frame), at a mean resolving power of 
700. The data allow us to measure several of the important Lick indices. The stellar populations and gas ionization 
state of the host galaxy of HE 1503+0228 are very similar to the ones measured for normal non-AGN galaxies. 
Dynamical information is also available for the gas and stellar components of the galaxy. Using deconvolution 
and a deprojection algorithm, velocity curves are derived for emission lines, from the center up to 4" away from 
the nucleus of the galaxy. Fitting a simple three-components mass model (point mass, spherical halo of dark 
matter, disk) to the position-velocity diagram, we infer a mass of M(r<lkpc) = (2.0 ± 0.3) 10 10 M within the 
central kiloparsec of the galaxy, and a mass integrated over 10 kpc of M(r<10kpc) = (1.9 ± 0.3)10" M Q , with 
an additional 10% error due to the uncertainty on the inclination of the galaxy. This, in combination with the 
analysis of the stellar populations indicates that the host galaxy of HE 1503+0228 is a normal spiral galaxy. 

Key words. Galaxies: dynamics, stellar populations - quasars: individual: HE 1503+0228 - techniques: deconvo- 
lution 



1. Introduction 

Luminous quasars are now known to be generally located 
at the cores of large galaxies. With the recent discovery 
that supermassive black holes appear to reside in most if 
no t all g alaxies with substantial bulges (e.g., Magorrian et 
al. 1998[ ), it is now quite likely that quasar-like activity is 
an extremely common but transient phenomenon, linking 
the growth of now extinct black holes to the violent pro- 
cesses during phases of active accretion (cf. McLeod et al. 
1999) ). However, very little of the physical processes oper- 
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* based on observations made with ANTU/UT1 at ESO- 
Paranal observatory in Chile (program 65.P-0361(A)), and 
with the ESO 3.5m NTT, at La Silla observatory (program 
62.P-0643(B)). 



ating during such episodes of nuclear activity is actually 
understood: neither do we know the actual conditions for 
fueling (or refueling) massive black holes, nor the time- 
scales involved. A particularly interesting problem is the 
importance of the feedback from a high-luminosity Active 
Galactic Nucleus (AGN) onto its host galaxy, in the form 
of huge quantities of ionizing radiation as well as possible 
mechanical outflows (jets). 

Quasar host galaxies have been studied almost ex- 
clusiv ely by imaging (e.g 



1997, Stockton 



Bahca ll et al. 

et al. |1998| , Marquez et al |200l|) . Recent Hubble Space 
Telescope optical studies have established that high- 
luminosity quasars generally reside in big ell iptica ls, irre- 
spective^ radio propertie s (Mc Leod & Rieke 1995 , Disney 
et al. |1995| , Hughes et al. |2000| ). There also seems to be a 
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Fig. 1. Left: field of view of 5'x 5 'around HE 1503+0228. The exposure (_ff-band) is 480 sec long. It has been obtained 
with SOFI on the ESO 3.5m NTT. The seeing is 0.8" sampled with a 0.29" pixel. Right: part of the VLT/FORS1 
pointing image. The seeing is 0.62"on this 30sec i?-band exposure. The 1" slits used to obtain the spectra of the target 
and of the PSF star are indicated. 



H 



R 



♦ 



♦ 



2 arcsec 



2 arcsec 




Fig. 2. Left and middle: deconvolved NTT and VLT images, where the quasar has been removed from the data. Right: 
zoom on the host galaxy and position of the slitlet used to obtain the spectroscopic data. 



trend that more luminous QSOs are hosted by more mas- 
sive galaxies (McLure et al. |l999j McLeod & Riekc |l995| ). 



A more detailed understanding of the physical condi- 
tions in host galaxies can only be obtained from spec- 
troscopic observations. Available evidence in this field 
is extremely scarce, and although extensive quasar host 
spectroscopy was conducted already in the early 1980s 
(e.g., Boroson et al. 1985), they have never really been 
followed up with improved instrumentation and analy- 
sis techniques, except for a few isolated objects (among 
them 3C 48 being the probably bes t-stu died case in 
this field 



eg 



Chatzichristou et al. 199£, Canalizo & 



Stockton |2000J) . Particularly important is the fact that 
nearly all these observations up to now were designed as 
"off-nuclear" spectroscopy, avoiding the strong contami- 



nation from the AGN itself, but yielding only information 
on the host for distances of <; 5 — lOkpc from the nucleus. 

With the aim of studying the stellar and gas content 
of quasar host galaxies, as well as their dynamics, we 
have initiated a spectroscopic campaign with the ESO- 
VLT to observe a sample of luminous radio-quiet quasars 
at low redshift. The sample is based on the Hamburg/ESO 
Survey (HES; Wisotzki et al. 2000 ), a wide-angle (~ 
7500 deg 2 ) search for optically bright (B ^ 17.5) quasars 
in the southern sky. This survey is particularly suited to 
provide samples for host galaxy studies, mainly because of 
two reasons: (1) because of its bright limiting magnitude, 
the survey yields low-redshift quasars in large numbers, 
(2) the selection does not encompass morphological crite- 
ria, i.e., there is no limitation to point sources as in most 
other optical quasar surveys. This fact and the wide range 
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Fig. 3. One dimensional deconvolved spectrum of HE 1503+0228 and its host galaxy. Each panel corresponds to one 
FORS1 grism and displays the individual flux calibrated spectrum of the quasar and its host alone. Note that the host 
spectrum shows no trace of the AGN broad emission lines, that are seen narrow in the host galaxy. Note also the very 
good agreement in the deconvolution of the 3 grisms in the overlapping regions. 



of selection criteria employed in the HES (cf. Wisotzki et 
al. 2000) ensure that quasars are selected irrespective of 
morphological properties of their host galaxies (see also 
Kohler et al. |1997| ). Our objects all have M B < -23 and 
z < 0.33. 

The present paper is aimed at exposing the techniques 
used to obtain high quality spectra of quasar hosts, decon- 
taminated from the contamination by the quasar, and to 
show their application to one object, taken as a test case: 
HE 1503+0228, at z = 0.135. 



2. Observations - data reduction 

2.1. Imaging 

We have obtained iJ-band images of HE 1503+0228, at 
La Silla observatory (ESO, Chile) on the night of 1999 
February 24 . SOFI was used, the near-IR imager of the 
3.5m ESO New Technology Telescope (NTT). The pixel 
size was 0.29" and numerous dithered images of the field 
were taken in order to produce a combined sky subtracted 
frame with a total exposure time of 480 sec and a final 
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mean seeing of 0.8". The image is displayed on the left 
panel of Fig. [l]. It shows the environment of the quasar on 
a 5'x 5' field, with no obvious companion. 

Additional data were obtained with the first of the four 
8.2m telescopes (UT1) at Paranal observatory. A shallow 
30s i?-band image was taken, in order to design the spec- 
troscopic mask. The seeing was 0.62", sampled with a pixel 
scale of 0.2". 

Both the near-IR and optical images were deconvolved 
using the MCS deconvolution algorithm (Magain et al. 
1998). While one single optical image of HE 1503+0228 
was obtained, many dithered IR images were available in 
the .ff-band. Eight stacks of images were produced with 
the data and simultaneously deconvolved as explained in 
Courbin et al. (1998). This procedure allows for better 
decomposition of the data into a sum of point source and 
"background" channels, hence a better extraction of the 
image of the host, uncontaminated by the AGN's light. 
The deconvolved data are shown in Fig. 0, where the 
quasar has been subtracted. 

The photometry of the host of HE 1503+0228 has been 
performed on the deconvolved images, in a 10" diameter 
aperture. The magnitudes we infer in that way are R = 
16.78 ± 0.20 and H = 14.57 ± 0.20. 

Determining the shape of the galaxy is of importance 
as the inclination of the disk component in the mass model 
(see next sections) depends on the observed galaxy shape, 
hence affecting the mass determination of the galaxy. We 
make the assumption that the galaxy has an intrinsic cir- 
cular geometry, to turn the fitted cllipticitics into an in- 
clination angle. Given the relatively low signal-to-noise of 
our images we chose to determine the shape of the host by 
fitting analytical profiles to the data (Kuhlbrodt et al. in 
preparation). Our model is composed of a nucleus, a disk, 
and a spheroid, convolved with the observed PSF before 
any fitting is performed. In the case of HE 1503+0228 the 
spheroid is used to model a bar-like feature in the center of 
the galaxy, that may affect the fit of the disk component 
if it were not taken properly into account. 

The analysis yields the Position Angle (PA=38° ± 3°) 
of the host, its inclination (i = 32° ±2°), and the angle 
(<p = 12 ± 3°) between the slit and the major axis of the 
galaxy. These values are in fact very close to the values 
found just by measuring the outmost isophotes in the VLT 
images. 

2.2. Spectroscopy 

The spectroscopic data were taken on the night of 2000, 
April II th with FORS1 mounted on Antu, the first of the 
four 8.2m VLT Unit Telescopes at Paranal observatory 
(ESO, Chile). Three grisms were used (G600B, G600R, 
G600I) to cover the whole optical spectral range 4000- 
9500A with a mean spectral resolution of about 700. The 
seeing during the observations was good, about 0.6". Each 
grism was exposed 1200 sec, under photometric conditions 
and dark time, which yields a mean S/N of about 150- 



200 per pixel (quasar+host). The large pixel size of the 
detector was used, 0.2". 

We show in the right panel of Fig. [I] and in Fig. || the 
orientation of the slitlets relative to the quasar host and to 
the PSF star. In Fig. |, the images have been deconvolved 
and the quasar has been removed from the data. The 
FORS mask consists of 19 slitlets with fixed length and 
width (19"x 1"). One slit was centered on HE 1503+0228, 
one slit was centered on a bright PSF star and the rest of 
the slits were used either to observe fainter PSF or galax- 
ies. Galaxy Gl is indicated in Fig. |l]. It is at z ~ 0.3 and 
is unrelated to the quasar. 

At the adopted spectral resolution, the area accessible 
on the CCD without loosing spectral information at the 
edges of the array is restricted to a strip of sky about 40" 
wide. As we have to observe simultaneously the quasar 
and at least one PSF star, the PA of the mask is there- 
fore imposed to us by the instrumental design, i.e., the 
40" strip of sky has to include both the quasar and the 
PSF star(s). This means that the slitlet centered on the 
quasar is not necessarily oriented along one of the axes 
of symmetry of our targets. It also means that we do not 
observe at parallactic angle. This is however not a major 
problem as FORS1 has an efficient atmospheric refraction 
corrector and the data remain unaffected by slit losses. 

A two-dimensional wavelength calibration was ap- 
plied in order to correct for slit curvature, followed 
by a two-dimensional sky subtraction. The spectrum of 
HE 1503+0228 and of the PSF star were reduced exactly 
in the same way, and they were rebinned to a common 
pixel size of lA in the spectral direction. They were also 
rebinned to a common starting wavelength, which ensures 
that to each spectral resolution element of the quasar, cor- 
responds the same spectral resolution element in the PSF. 

3. Extraction of the host's spectrum 

Most current spectroscopic studies of quasar host galaxies 
are designed as off-axis observations, with the slit of the 
spectrograph placed one to several times the size of the 
seein g disk away from the center of the AGN (e.g., Hughes 
et al. |2000| ). 

While such an observational strategy is minimizing 
contamination by the quasar, it has severe drawbacks. 
First, one can only study the outer parts of the galaxy. 
This is limiting the study to particular regions of the 
hosts. It also requires very long exposure times on large 
telescopes to detect the extremely faint wings of the hosts. 
Second, all velocity information is lost, as is the spatial in- 
formation in general, along the galaxy's axes of symmetry. 
Finally, the method is still restricted to the most extended 
objects and does not allow a full decontamination of the 
quasar's light. 

Instead, our observations are taken on-axis, placing the 
slit of the spectrograph on the quasar's center of light. 
This avoids all drawbacks of the off-axis observations, but 
requires clean post-processing techniques in order to re- 
move accurately the spectrum of the AGN. For this pur- 
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pose we use the spectroscopic version of the MCS deconvo- 
lution algorithm (Magain et al. 1998, Courbin et al. 2000| ). 
Using the spatial information contained in the spectrum 
of PSF stars located close to the object of interest, the 
algorithm sharpens the data in the spatial direction, sub- 
samples them in order to achieve smaller pixel size, and de- 
composes them into two channels containing (i) the spec- 
trum of the AGN and (ii) the two-dimensional spectrum 
of the host galaxy. 

In the deconvolved spectrum, the new sub-sampled 
pixel size is 0.1" and lA in the spectral direction. In Fig. || 
are displayed the one dimensional flux calibrated spectra, 
integrated along the spatial direction on the deconvolved 
data. Interestingly, the three FORS1 grisms have signifi- 
cant regions of overlap that are in perfect agreement with 
one another. In addition, at the redshift of HE 1503+0228, 
most emission lines are present in two grisms. This allowed 
us to construct a combined 2-D flux calibrated spectrum. 
All measurements presented in this paper, except for ve- 
locity curves, are done on the final combined spectrum. 



4. Dynamics of the host galaxy 

4.1. Redshift 

From the different emission lines of the quasar and of the 
galaxy, we can extract a sharp value for the redshift of 
HE 1503+0228, by simple Gaussian fitting. The lines of 
the host galaxy in the deconvolved spectrum have to be 
measured with particular care: in the center of the host, 
the lines are double, due to the velocity field of the galaxy 
and to the deconvolution method. This is mainly due to 
the extreme contrast between the narrow emission lines of 
the host, and the broad lines of the quasar: some of the flux 
in the sharpest parts of the host's emission lines is taken 
in the point-source component (the quasar) of the decon- 
volved spectrum, i.e., the central parts of the lines are 
dimmed. We have therefore taken the mean wavelength 
between the two sides of emission lines, when measuring 
the host's redshift. The results for each line are given in 
table [|. The details of the individual emission lines is ac- 



Table 1. Redshift of HE 1503+0228, as infered from the 
emission lines of the host galaxy alone, and from the 
quasar alone. The redshifts are the same for the quasar 
and for its host. 



cessible via a method described in Section 4.2, which does 
not alter the center of the lines. 

4.2. Velocity curves: extraction method 

In the spectroscopic version of the MCS algorithm 
(Courbin et al. |2000| ), the difference in spatial properties 
between the AGN (point source) and the host galaxy (ex- 
tended component) is used to separate the spectra of these 
two sources. A smoothing in the spatial direction is applied 
for this purpose to the extended component. While such 
a decomposition method is well suited to most spectral 
regions of the host, it is not optimal for deriving velocity 
curves: spatial smoothing (which is not precisely known as 
it varies with the local S/N) may modify the spectral po- 
sition of the lines at a given spatial position, by averaging 
spatial components of different radial velocities. 



Line 


rteasmrt 


rteasmrt 




(host only) 


(quasar only) 


Oil (3727A) 


0.1353 


0.1352 


H/3 (486lA) 


0.1357 


0.1352 


OIII (4959A) 


0.1362 


0.1349 


OIII (5007A) 


0.1356 


0.1348 


OI (6300A) 


0.1355 


0.1352 


Nil (6548A) 


0.1352 




Ha (6563A) 


0.1353 


0.1352 


Nil (6583A) 


0.1354 




SII(6716A) 


0.1356 


0.1343 


Mean 


<0.1355> 


<0.1350> 


Mean error 


(±0.0003) 


(±0.0003) 



We have devised another method which is better suited 
for extracting radial velocity curves of the host galaxy. 
Instead of using the difference in spatial properties to ex- 
tract the spectra, we take advantage of the different spec- 
tral properties of the AGN and host to determine the 
radial velocity as a function of distance. 

The fact that the spectral lines of the host are much 
narrower than any feature in the AGN spectrum is used 
to determine the spectral position and intensities of the 
host spectral lines as a function of spatial position. More 
specifically, at each spatial coordinate we fit the sum of 
a Gaussian profile of variable width representing the host 
emission line, plus a smooth numerical curve accounting 
for the AGN spectrum plus the continuum of the host. 

In doing this fit, we also take advantage of the con- 
stant shift in wavelength between two given emission lines: 
several emission lines are fitted simultaneously under the 
constraint that they are separated by a fixed gap in wave- 
length. Thus, when we refer to the Ha velocity curve, we 
in fact use the Ha and [Nil ] emission lines simultaneously. 
Similarly, the H/3 curve is obtained from the H/3 and the 
[OIII] (5007A) lines, simultaneously . An example of such 
an extraction is displayed on Fig. ^. 

In converting the wavelength shifts into velocity 
curves, we take several effects into account: 

(1) The orientation of the slit does not coincide with the 
major axis of the galaxy. This is due to observational 
constraints, i.e. the fact that we had to obtain the spec- 
trum of a PSF star simultaneously with the quasar and 
host (see Section 2.2). 

(2) The slit width is comparable to the angular size of the 
galaxy. The observed wavelength shifts are thus aver- 
ages of velocity components coming from many parts 
of the galaxy. 

(3) The galaxy is not seen edge-on, but with an inclination 
angle i. Radial velocities have thus to be divided by 
sini to be transformed into rotation velocities. 
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Fig. 4. Example of emission lines extraction. Here are the 
Ha and [Nil ] emission lines in /-band grism. The left part 
shows a zoom on initial reduced spectrum. The middle 
one shows the emission lines alone whereas the right part 
shows the quasar and galactic continuum. 

(4) The radial velocity curves are smeared by convolution 
with the seeing profile, both in spatial and spectral 
directions. 

The procedure we use assumes a simple analytical 
mass model for the host galaxy, which includes a thin 
disk, a spherical dark matter halo, and a central point 
mass. The velocity at a distance r away from the galaxy's 
center is then given by: 



V mo d{r) = \/-r(Fdi S k + Fdark + F C m) 



(1) 



The radial forces in equation 
a Kuzmin potential for the disk: 



-Fdisk(r, z) 



r 2 GM disk 



(r 2 + (a + \z\) 2 Y- 



111) are computed using 



(2) 



where we assume z = 0, hence neglecting the thickness of 
the disk. Mdisk is the total mass of the disk and a measures 
its scale length. 

Dark matter is assumed to be spherically symmetrical 
about the galaxy's center. The force it creates at the dis- 
tance r away from the center, given a radially dependent 
density distribution is: 



-Fdark(r) 



( arctan — ) 

ro T-o 



(3) 



where ro is the scale length for the mass distribution 
within the dark matter halo and po its central density 
(Binney & Tremaine [l987|) . 

The force from the central mass M cm is simply 



Fcu(r) = - 



GAL 



(4) 



The five parameters Mdisk, a, Po, r$ and M cm are de- 
termined by a least squares fit to the data, i.e., to the 
measured velocity shifts at each spatial pixel along the 
slit. We therefore minimize the following % 2 : 



N 1 



^riican }j 



(5) 




Major axis 




Slit axis 



Fig. 5. Schematic illustration of geometrical effects con- 
sidered for deriving the velocity curves. Top : the original 
face-on galaxy is assumed to be circular (dotted line) . Seen 
with an inclination i, it appears elliptical (dashed line). In 
addition, the slit is not aligned with the major axis of the 
inclined galaxy. Bottom : The mean velocity at one po- 
sition along the slit, which would be integrated along a 
line perpendicular to the slit (dashed-dotted line) if the 
galaxy was seen face-on, has now to be integrated along 
a line which is inclined (and lengthened) with respect to 
the slit (dotted line; see text). 

where N is the total number of pixels along the slit, and 
where dj ± o~j is the measured radial velocity shift at pixel 
j along the slit. F mean is a one dimensional vector that 
contains the mean (see next paragraph for the significa- 
tion of "mean") velocity predicted by our mass model at 
pixel j along the slit. Due to the seeing, this velocity is af- 
fected by blurring in the spatial direction, hence the need 
for the convolution by the 1 dimensional PSF s obtained 
from the spectrum of stars. As the galaxy has an inclina- 
tion i, we multiply the result of the convolution by sin i 
before comparing with the observed velocity dj . The error 
o-j on the radial velocity measurements is estimated by 
comparing the measured dj for the same emission lines in 
different grisms. This was possible for the Ha, [Nil], H/3 
and [OIII] emission lines. 

All the calculations in Eq. || assume that the galaxy is 
seen face-on, the inclination being taken into account by 
the sinz term. Due to the non-zero angle <fi between the 
slit and the major axis of the galaxy, the velocity in Eq. || 
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Fig. 6. Description of the integration domain for the com- 
putation of V mcan (P). The major axis of the galaxy is 
along X. / e and fi are the limits of the integral also 
shown in Fig. Only the component of the model ve- 
locity Vmod along the Y axis contributes to the integral 
when the galaxy is seen inclined. 

is a weighted mean V 1Ylcan of several velocity components 
across the slit. 

For each pixel j along the slit, lets say at point P (see 
Fig. H), one can compute a mean velocity as follows: 

Kncan(P) = [ I[r(x)] ■ V nlod [r( X )} ■ ■ dx (6) 

hot J X{ . r{x) 

V mean (P) is the average of all the velocity components 
located along the dash-dotted line in Fig. If the galaxy 
was seen face-on (as we would like it to be in Eq. [3|), this 
line would be rotated and lengthened, and mapped onto 
the dotted line of Fig. ||. The integral in Eq. ^ runs along 
this dotted line. Its equation (Eq. 0) is easily obtained 
from geometrical considerations, given the inclination i of 
the galaxy with respect to the plane of the sky and given 
the orientation <j> of the slit relative to the major axis of 
the galaxy: 

coscj) x P 

y = — - (x -)• (7) 

cos I cos <p 

xp is the coordinate of the point P on the major axis 
X (Fig. ^J) of the galaxy. 

The integration in Eq. ^ is then performed as follow: 
for each position x along the major axis of the galaxy, 
we calculate the distance r{x) from the galaxy's center 
O to the point of the integration line that is considered. 
This distance is used to weight each velocity component 
Knod[^(a;)] of the velocity model, according to the light 
distribution of the galaxy in the emission line. We have 
estimated the light distribution I(r) of the galaxy, assum- 
ing it is intrinsically symmetrical about the center of the 
galaxy. We model it by fitting a Moffat profile to the 1- 
dimcnsional intensity profile along the emission line, and 
applying deprojection (note that the fit of the emission 
line is one-dimensional, but since I(r) is known analyti- 
cally, one can compute weights in 2 dimensions when cal- 
culating the integral in Eq. fy. When doing the fit, the 



400 



-400 





1 1 1 


1 ' 










v(radius) 




Ha 






- 


Grism I 


I , 




- 






• 










- 




• 

I 


- 






| 










■ data 








! 3 1 






, 1 , 





-lxlO 4 -5000 5000 10 1 

Radius (pc) 



Fig. 7. Inclination/seeing corrected velocity curves for the 
Ha and [Nil] emission line (/-grism only): the dots show 
the velocity measurements. The solid line superposed is 
the fit of our three components model after taking into 
account misalignment between the slit and the major axis 
of the galaxy (see text) . This solid line shows the variation 
of Vmcan along the slit. It is then corrected for inclination 
and seeing, following Eq. M, to give the true velocity as a 
function of distance to the center of the galaxy. 



Moffat profile is convolved with the PSF. Correct nor- 
malisation is assured by /tot : the total flux of the regions 
contributing to the measured emission at position P along 
the slit. 

The term in Eq. ^| is the projection of the rota- 
tion velocity on the minor axis of the galaxy, which will 
give rise to a radial velocity once inclination is taken into 
account. The component of the velocity along the X-axis 
does not contribute to the mean velocity calculation, once 
inclination has been taken into account. 

The intersections between the integration line and the 
slit edges, fi and f e , give us the limits of the integration, 
as shown in Figs. || and ^). 

An example of a inclination/seeing-corrected velocity 
curve is given in Fig. |?] for the Ha/N[II] lines. In our 
calculations of the linear scales and distances we have used 
Z/o=65 km.s _1 Mpc _1 , A = 0, and Qm = 0.3, leading to a 
scale of 1" = 2.5 kpc. 

4.3. Velocity curves: results 

We give in Table || the best fit for the 5 parameters of our 
mass model. The mean velocity curve, as extracted from 
the Ha and H/3 lines is displayed in Fig. ||. This mean 
corresponds to four independent curves obtained for the 
Ha and the H/3 lines, both seen in two different grisms. 
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Fig. 8. Mean velocity curve between four models fitted 
on the Ha line (2 grisms) and on the H/3 line (2 grisms). 
The error bars are derived from the dispersion between 
the four curves. The result is corrected for seeing and slit 
effects, as well as for inclination. 

Table 2. Mass model parameters for our velocity curves 
in each individual grism and for each emission line, a and 
r are in parsec, po in M /pc 3 and Mdisk and M cm in M . 



Line 


Grism 


a 


Mdisk 


Po 


ro 


Mom 


Ha 


R 


1088 


7.9 10 iu 


0.30 


2033 


6.1 10' 


Ha 


I 


1163 


8.6 10 10 


0.31 


2074 


4.8 10 7 


H/3 


B 


908 


4.8 10 10 


0.56 


1753 


3.5 10 6 


H/3 


R 


1428 


5.9 10 10 


0.23 


2196 


2.5 10 5 



Table 3. Weighted means of the masses obtained in Table 
U where the weights are computed from the residuals of 
each individual fit. M to t is the sum of Mdisk and Mdark- 
Units are solar masses. 



Masses 


r < 1 kpc 


r < 3 kpc 


r < 10 kpc 


M disk 


1.9±0.3 10 IU 


5.1±0.7 10 1U 


7.0±0.9 10 1U 


Mdark 


1.1±0.3 10 9 


1.6±0.3 10 10 


1.2±0.1 10 11 


Mtot 


2.0±0.3 10 10 


6.7±1.0 10 10 


1.9±0.2 10 11 



Table || gives our mass estimates which are weighted 
means of the four fits presented in Table ||, the weights 
being a function of the residuals after modeling of the 
data obtained in each individual grism. These weights are 
the same as the ones applied to compute the mean velocity 
curve of Fig. [|. 

We investigate in Table [| the effect of the error on the 
inclination of the galaxy and on its PA. Uncertainties on 
inclination and orientation of the major axis lead to error 
bars on masses, which are summarized in Table m. 



Table 4. Effect of the error on the inclination (32 ± 2°) 
and PA (38 ±3°) of the host galaxy. Results are given for 
weighted means of 4 velocity curves, as in Table |, varying 
the inclination and the PA within their error bars. 



Masses (M©) 


r < 1 kpc 


r < 3 kpc 


r < 10 kpc 


Inclination 
Mdisk 

Mdark 


1-9+0.2 10l ° 


5.1;™ io 10 
i.6jo°:S io 10 


7.0-11 10 10 
1.2;g;J 10 11 


M tot 


2.0;^ io 10 


6.7;- 8 7 10 10 


i-9;L 2 io 11 


Major axis 
Mdisk 

Mdark 


1 9+ - 02 10 10 

1 1+ 01 10 9 
1 - 1 +0.04 lu 


5.1±g'i 10 10 

i.6±g;J a io 10 


7.0+° 2 2 10 10 
1 2+ 001 10 11 


Mtot 


2 0+ 03 1 10 




1 q+0 03 10 n 



These results should be interpreted with care, as there 
is a strong correlation between the model parameters and 
since the adopted model is simpler than reality. However, 
while the details of the mass distribution within the galaxy 
remain inaccurate, we can still determine the total mass 
of dark matter and as well the disk mass enclosed in a 
given radius. 

Given the results in Tables ||, || and ||, we infer a mass 
estimate of M(r < 1 kpc) = (2.0±0.3)10 10 M . The corre- 
sponding mass, integrated over 10 kpc is M(r < 10 kpc) = 
(1.9±0.3)10 n M . Table I shows that an additional error 
of 10% should be added, due in part to the uncertainty on 
the PA of the major axis of the host, but mainly to the 
uncertainty on its inclination. 



5. Stellar and Gas Content 

One of the primary goals of our study is to infer the char- 
acteristics of the stellar population of quasar host galaxies: 
does the stellar population of HE 1503+0228 match that 
of a normal galaxy ? Where does it lie in the Hubble se- 
quence ? 

Let us first note that H = 14.57 mag corresponds, with 
the same cosmological parameters used earlier, to an ab- 
solute magnitude of M# =—24.5 mag. This, in addition to 
the profile fitting described in the previous sections, indi- 
cates a normal early type local spiral galaxy (Bothun et al. 
1985; Arimoto & Jablonka 1991). While the photometric 
data alone already allow to point out a plausible morpho- 
logical type for the host of HE 1503+0228, spectroscopy 
is the key to the determination of its stellar content. A 
few considerations have to be made before looking at the 
spectral features in any detail. 

First, galactic extinction has to be taken into account. 
The extinction towards HE 1503+0228 is A B = 0.217 
corresponding to E(B-V)=0.05 (Schlegel et al 1998). Our 
spectrum of HE 1503+0228 has been corrected accord- 
ingly. 

Second, although no velocity standard stars were taken 
during the observations, the resolution in the i?-grism 
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Fig. 9. The relation between the equivalent widths of 
Ha+[NII] and [Oil]. The black dots show the normal 
spiral galaxies of Kennicutt's sample. The plain line cor- 
respond to [Oil] = 0.4 Ha+[NII], mean relation of the 
galaxies dominated by photoionization. The cross indi- 
cates the position of HE 1503+0228. 

(4.9A; 3.6A at rest-frame) is good enough to make use 
of the stellar library of Jacoby et al. (1984) (FWHM =3.5 
A). This gives a velocity dispersion of ~ 170 km/s, which 
we only take as an indicative upper limit value. 

Third, we observed HE 1503+0228 with three different 
grisms, each having a different dispersion. The three spec- 
tra were combined in order to cover the wide wavelength 
range 3670A to 9165A (3150A to 8073A once corrected 
from the redshift z=0.1355, or v=37837.4 km/s). In doing 
this, we have degraded the observations to the resolution 
of the less resolving grism, i.e., 5.7 A (or ~ 5 A at rest- 
frame) . 

Finally, due to the low signal-to-noise ratio of the spec- 
trum in its bluer part, we consider only features from 4000 
A on (3523A at rest-frame). 

We have only considered the integrated spectrum of 
HE 1503+0228. In fact, the signal-to-noise of the mea- 
surements on each individual spatial resolution element of 
the 2D data varies between 3 and 6, and therefore does 
not give much hope in tracing any radial variation in the 
stellar population, unless deeper data are obtained. 

In the following, we measure well known indicators 
for the underlying stellar populations, for the host of 
HE 1503+0228, and we compare them with the values 
obtained for two control samples: the one of Trager et al. 
(1998) for the elliptical galaxies, and the one of Kennicutt 
(1992a, b) for the spiral galaxies. 

5.1. Emission lines 

One of our concern was to verify that the deconvolution 
procedure had indeed left out any non thermal excitation 
process. 



Fig. 10. The relation between the 4000A break and 
the equivalent width of H/3 for normal spiral galaxies 
(Kennicutt 1992) (dots) and HE 1503+0228 (cross). 



We have measured the equivalent widths (EW) of the 
most prominent emission lines seen in the spectrum of 
HE 1503+0228, i.e., [Oil], [OIII] (4959A), [OIII](5007A) 
andHa+[NII] (see Table |). 

Kennicutt (1992a) provided measurement of the 
Oxygen lines [Oil] (3727A), [OIII] (5007A) and 
Ha+[NII]. He noted that among galaxies with nor- 
mal star formation rates (EW(Ha+[NII]) < 40A), the 
[OIII] lines are rarely detected at all (EW < 0.5A). For 
HE 1503+0228, while EW(Ha+[NII]) - 30A, we detect 
[OIII], with an equivalent width of ~ 4A. This is the only 
and discreet sign of an excitation process in or surrounding 
HE 1503+0228 which is not directly related to star for- 
mation and, hence, which might be triggered by the cen- 
tral AGN radiation field. Kennicutt (1992a) also reported 
that galaxies dominated by stellar photoionization follow 
a mean relation EW([OII]) - 0.4 EW(Ha+ [Nil]) (with 
an rms of 50% ). HE 1503+0228 falls nearly exactly on 
this relation, as shown in Fig. |[ Kennicutt's galaxies are 
represented by plain dots, the mean relationship between 
Ha+[NII] and [Oil] is the solid line. The position of 
HE 1503+0228 in this plane is indicated by a cross. 

HE 1503+0228 exhibits a weak [OIII] (5007A) line. 
We have left out from the Kennicutt sample the ob- 
jects which displays strong [OIII] (5007A) lines, when 
analysing in the next section the absorption line features. 
This ensures that we compare HE 1503+0228 with objects 
of similar properties. 

To further investigate this matter, we have measured 
the line intensity ratios, [OIII]/ H/3 (0.05), [OKI]/ H/3 
(0.54), [Nil]/ Ha (-0.20). Due to the deconvolution proce- 
dure, the ratios between Oxygen and Balmcr lines are up- 
per limits. Indeed, the Balmer lines are probably underes- 
timated. When placed in the diagnostic diagrams of Sodre 
and Stasihska (1999), the host galaxy of HE 1503+0228 
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Table 5. Lick indices for the spectral features in the host galaxy of HE 1503+0228. 



Name D4000 G4300 Ca4455 Mg 2 Fe5270 Fe5335 [QUI] Hq+ [Nil] 
Value L48 3A1 2~M 0.157 L89 L12 VI 30 
Error 0.04 0.16 0.12 0.002 0.08 0.11 0.3 2. 




Fe5270 

Fig. 11. The relation between the Fe5270 and Mg 2 index 
for the normal spirals of Kennicutt et al. (1992) (dots) 
and the ellipticals of Trager et al. (1998) (open triangles). 
As in Fig. | and [l(], the cross indicates the position of 
HE 1503+0228. 



falls perfectly on the sequence drawn by the normal spiral 
galaxies. 

To summarize, our analysis of the emission lines sug- 
gests that HE 1503+0228 is a perfectly normal spiral 
galaxy, where the gas emission is dominated by ionization 
by stars rather than by the central AGN. 

5.2. Absorption lines 

Our spectra have sufficient signal-to-noise to measure the 
Lick indices for several stellar absorption features. As the 
spectral resolution of the Lick system is lower than the 
resolution of our data, we have convolved the spectrum of 
HE 1503+0228 with the appropriate Gaussian kernel, in 
order to match a 200 km s _1 rest-frame resolution. 

For a full compatibility between the measurements 
done for HE 1503+0228 and for the galaxy control sam- 
ple (see previous section), we have re-measured the Lick 
indices on Kennicutt 's galaxies in exactly the same way 
as for HE 1503+0228. 

A difficulty when looking at the stellar population in 
HE 1503+0228 is to avoid contamination due to the gas 
emission. This emission can either "fill-in" absorption lines 
-the Balmer lines is a well known example, but Mgb can 
be affected as well (Goudfrooij & Emsellem 1996)- or con- 
taminate the continuum side bands used for the measure- 
ments of the equivalent widths. We have identified four 
absorption lines that are strong enough and free of this 



contamination by emission lines. They correspond to the 
Lick indices G4300, Ca4455, Ca4668, Fe5270 and Fe5335. 
With only a tiny change in the definition of the contin- 
uum window used for the Mg2 index, we could avoid con- 
tamination by the [OIII] (4959A) emission line for this 
index. The small change we impose on the continuum side 
bands produces a change of at most ±0.04mag (or 20%) 
in the measurement of the Mg2 index. In addition to the 
indices for the absorption lines, we also measure the im- 
portant calcium break at 4000A, which gives the D4000 
index (Bruzual 1983). 

Table | summarizes the results of our measurements. 
The available signal-to-noise ratio (per A) for performing 
the measurements varies between ~ 10 (blue part) and 
~25 (red part). 

The continuum shape of HE 1503+0228 is rather flat, 
indicative of on-going star formation. Similarly to the re- 
sult of Fig. g Fig. [ft] puts HE 1503+0228 among nor- 
mal spiral galaxies. The latter figure compares the index 
D4000, which measures the 4000A break, with the equiv- 
alent width of H/3 which is indicative of the presence of 
star formation. HE 1503+0228 nicely falls on the sharp 
relation traced by normal spirals. 

As for the chemical abundances, we chose to present 
in Fig. |l| the Mg 2 and Fe5270 indices. They trace 
the a and iron peak-elements, respectively. Here again, 
HE 1503+0228 is located exactly on the Hubble sequence, 
at the location of spiral galaxies, rather than ellipticals. 
Furthermore, the values of the absorption line equivalent 
widths strongly suggests the stellar population of an early 
type spiral galaxy. 

6. Conclusions 

We have undertaken a VLT program aimed at unveiling 
the spectroscopic properties of quasar host galaxies for 
a sample of radio quiet, bright quasars at low redshift. 
Our long term scientific goal is double: (1) to compare 
the stellar populations of quasar hosts with those of other 
("normal") galaxies, and (2) to study their dynamics, as 
close as possible to the central AGN. 

The present paper considers HE 1503+0228 as a test 
case. An important part of our work is dedicated to the 
correction of geometrical effects and removal of the atmo- 
spheric blurring. This is crucial in the central parts of the 
galaxy, with sizes comparable to the seeing disk. 

We find that the stellar population of the host galaxy 
of HE 1503+0228 compares well with that of the nor- 
mal non- AGN spiral galaxies of Kennicutt (1992a, b) and 
Trager et al. (199S). They show no trace of enhanced star 
formation. The interstellar medium does not show sig- 
nificant ionization by the central AGN. Its mass M(r < 
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10 kpc) = (1.9 ± 0.3)10 n M Q also compares very well 
with that of other spiral galaxies. To summarize, the host 
galaxy of HE 1503+0228 is a normal spiral galaxy. 
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